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0^ ' Abstract. - The study of the ground state of spinless fermions in 2d disordered clusters 

(Phys. Rev. Lett. 83, 1826 (1999)) has suggested the existence of a new quantum phase 
for intermediate Coulomb energy to kinetic energy ratios r s . Exact diagonalization of the same 
small clusters show that its low energy excitations (quantum ergodicity above a few "hexatic" 
excitations characterized by oriented currents) significantly differ from those occuring in the 
Fermi glass (weak r s ) and in the pinned Wigner crystal (large r s ). The "hexatic" excitations 

i-rt ■ vanish for temperatures of order of the Fermi temperature. 
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An insulator-metal transition (IMT) has been observed when the carrier density n s is 
increased in dilute two dimensional (2d) gases of electrons [[j], || and holes || having in common 
very low Fermi temperatures (Tp «1 K). The critical Coulomb energy to Fermi energy ratios 
r s oc 1/ ' ^JrTs are typically around 10 — 35 at the IMT, dropping as a function of the elastic 
scattering rate 1/r from r s rs 35 for a very clean heterostructure |3|] to a value of order 10 when 
r < 10 _11 s. This coincides with the critical values rf found in numerical simulations for the 
melting of the Wigner crystal, in the clean limit (r^ rs 37) and in the disordered case (r^ ~ 10, 
see Refs. || 0). In the insulating phase (r s > r^), I ~ V characteristics [|| || consistent 
with the existence of a pinned crystal have been observed, accompanied by a sharp drop of 
the compressibility [pi. Therefore, it has been suggested || [|, || that the insulating character 
of the phase is due to the pinning of a Wigner crystal, and not to Anderson localization. 
At higher densities (r s < rf rs 6), the kinetic terms become more important, giving rise 
to a re-entry where the metallic behavior disappears. In this limit, the insulating behavior 
predicted by the scaling theory of localization is recovered M (Fermi glass). The metallic 
behavior characterizing intermediate ratios rf < r s < r^ persists 0] down to T — > 35 mK 
(T/Tp — > 10~ 2 ), but vanishes H, B| when the temperature T exceeds the Fermi temperature 
Tp. Eventually, the IMT occurs when the elastic mean free path I and the Fermi wave vector 
kp are such that kpl rs 1. In summary, the interaction, disorder and kinetic effects are a priori 
relevant and non pertubative in those dilute gases and the observation of the IMT in different 
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materials leads us to suspect the existence of a rather general mechanism. 

Although a generalization up to r s — > 10 of the non-interacting theory where the disorder 
depends on temperature has been proposed to partly explain H the metallic behavior, we 
numerically explore another alternative: Would the IMT be one of the signatures of a new 
quantum phase (Coulomb metal) driven by the interplay between kinetic energy and Coulomb 
repulsion in a random substrate? A first numerical study H was focused on the ground state 
of spinless fermions with Coulomb repulsion in disordered clusters forming a 2d torus enclosing 
an Aharonov-Bohm flux <p along the longitudinal direction. The topology of the pattern of 
the driven currents provides the signature of a new intermediate quantum regime. At weak 
r s (Fermi glass) one has a 2d pattern where the local current flows in every direction. At 
intermediate r s (new regime), the pattern becomes essentially Id (no transverse currents and 
enhanced Id longitudinal flows enclosing the flux, see also Ref.|l). At large r s (pinned Wigner 
crystal), the persisting Id flows disappear and charge crystallization sets in. This suggests the 
possible existence of a new quantum phase characterized by an orientational order of the 
persistent currents without charge crystallization. This phase could be a quantum analog in 
a disordered system of the hexatic phase supposed to occur at the melting of the classical 



Wigner crystal 1 10 above the Kosterlitz-Thouless temperature. Characteristic values rf and 
Tj delimiting the new intermediate regime were obtained in agreement with the observed 
re-entry {rf « 4) and IMT {rf « 10). 

In this Letter, we study the low energy excitations of the same clusters which were studied in 
Ref. [0 , and their statistics when the microscopic configurations of the random substrate are 
changed. We find that the characteristic low energy excitations of the new intermediate regime 
disappear for excitation energies e of order of the Fermi energy ef- Inside the intermediate 
regime, the low energy excitations (e < ef) are characterized by a non random orientation of 
the local currents and deviations from Wigner-Dyson (W-D) spectral statistics. In analogy 
with the classical Coulomb problem, we refer to these oriented non chaotic excitations as 
"hexatic". For excitation energies e > ep, the local currents become randomly oriented and the 
levels obey W-D statistics (quantum ergodicity). In the Fermi glass and in the Wigner crystal, 
ergodicity is expected for energies much larger than ep and the levels remain uncorrelated 
for e rs ep. Quantum ergodicity with correlated energy levels for e > ef is then another 
characteristic of the intermediate regime, in addition to "hexatic" excitations for e < ef- 
Moreover, the sensitivity of the levels under a change of the threaded flux has a remarkable 
temperature dependence. When T — ► Tp, the enhancement observed in the intermediate 
regime is suppressed while the weak responses of the Fermi and Wigner limits are enhanced. 
Eventually, assuming that this intermediate regime observed in small clusters is the signature 
of a new quantum phase, we propose the local current j driven by a flux as a possible order 
parameter for the successive quantum transitions. The first transition would be associated to 
the angle of j with the direction enclosing the flux and the second one to its amplitude j. 

We consider a disordered square lattice with N — 4 spinless fermions on L 2 = 36 sites. The 
Hamiltonian reads 

ff—'E^ + E^ + ^ESr 1 ' (!) 

where c\ (c$) creates (destroys) an electron on site i = (i x ,i y ), the hopping term t between 
nearest neighbours characterizes the kinetic energy, the site potentials Vi are taken from a box 
distribution of width W, ri{ = c\ci is the occupation number at the site i and U measures the 
strength of the Coulomb repulsion, r-y is the inter-particle distance for a 2d torus. In our units, 
the factor r s = U(2ty/im e )~ l at a filling factor n e — N/L 2 . The Fermi energy is defined by 
ef = Aim e t and a Fermi golden rule approximation gives kpl ~ lQ2im e (t/W) 2 . We consider 



G. BENENTI et al.: Intermediate low temperature Coulomb 2d phase O 

W/t = 5 corresponding to kpl — 2.7. The boundary conditions (BCs) are always taken 
periodic in the transverse y-direction, and such that the system encloses an Aharonov-Bohm 
flux <f> = it/2 in the longitudinal x-direction (0 = tt corresponds to anti-periodic BCs). The 
Hamiltonian is diagonalized using the Lanczos algorithm for a statistical ensemble of 10 3 
clusters. The first 20 energy levels E n (n = 0, 1, 2, . . .) are considered for values of r s chosen 
inside the three phases (r s — 0.8 for the Fermi glass, r s — 6.3 for the conjectured intermediate 
phase and r s = 30 for the pinned Wigner crystal). 

We begin the presentation of our results by discussing the spectral fluctuations. For the 
one body spectra, localized wavefunctions yield uncorrelated spectra with Poisson statis- 
tics, characterized by a distribution P(s) of the energy spacings between successive lev- 
els going to Pp{s) — exp(— s) when L ^> L\. Delocalized wavefunctions yield correlated 
spectra with W-D statistics. With time-reversal symmetry (</> = 0), P(s) — > P\y( s ) = 
(irs/2) exp(— its 2 /A) as L increases, and without time-reversal symmetry (0 = 7r/2), P(s) — > 
Pwi s ) = (32s 2 /7r 2 ) exp(— 4s 2 /tt). The statistical properties of many body spectra are more 
complex. Without interaction, the levels are uncorrelated, and a sufficient two body interaction 
may restore W-D statistics. This is at the basis of the success of random matrix theory to 
describe complex nuclei, where low energy quasi-particle or collective excitations are followed 
by chaotic high energy excitations. The low energy modes being close to integrability display 
Poisson statistics as non-generic integrable models, whereas in the absence of integrals of 
motion quantum ergodicity and W-D statistics set in [ pjj . The question is to know at what 
energy threshold the Poisson- Wigner crossover takes place. An answer was proposed in Ref. 
[ [La H3] based on the breakdown of a perturbation theory valid for weak U (the matrix 
element of a Slater determinant to the "first generation" of determinants directly coupled 
to it by interaction is of order of the level spacing of the latter determinants). This question 
was recently addressed H, [L5| for the Hamiltonian which we consider, showing as expected 
quantum ergocity at sufficient excitation energies, but failing to check the proposed criterion 
(the estimation of the matrix elements deserving a particular study) . Another question is to 
know the factors r s delimiting the regime of quantum ergodicity for a given excitation energy 
e. The qualitative answer can be guessed from a simpler problem Jig , |l7| (two interacting 
particles in lei) where a duality t/U «-> U/t implies that a maximum spectral rigidity occurs 
for U s=s t. This leads us to a general picture where Fermi glass and Wigner crystal should 
be associated to integrability and Poisson statistics, in contrast to an intermediate regime 
displaying quantum ergodicity at low excitation energies. 

We study the distribution of the normalized TV-body energy spacings s n — (E n+ i — E n )/ < 
E n+ \ — E n > (the brackets denote ensemble average). In Fig. [j] A, one can see that P(s) is 
given by the Wigner surmises, with an orthogonal-unitary crossover when <j> = — ► 7r/2, in 
the intermediate regime for 1.4 < e/ep < 1.9 (n = 10, . . . , 20). We define a parameter 

= var(P( S ))-var(P^( S ))) 
1 var(P P ( S ))-var(P£( S )))' [) 

(var(P(s)) denoting the variance of P(s)). rj — 1 when P(s) = Pp(s) and r\ = when 
P(s) = P^r(s). Taking (f> = n/2, the variation of rj as a function of e/ep is given in Fig. p] B 
in the three phases. For the Fermi glass with r s — 0.8, the two first spacings display r/ values 
compatible with the unitary ensemble, while the following spacings become more and more 
Poissonian. This is what is expected for r s = 0, the two first excitations being single-electron 
transitions with W-D statistics for L < L\, while the following excitations being the sum 
of more than one single-electron excitation, give uncorrelated levels. For a small r s , this still 
remains true, the system being composed of weakly interacting quasi-particles. For the Wigner 
crystal [r s — 30), ij does not depend on e up to 3ep and keeps an intermediate value, rj = 1 
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Fig. 1. - A: Spacing distribution P(s) for r a — 6.3 when <j> — (empty circles) and <j> — 7r/2 
(full circles) at excitation energies 1.4 < e/eF < 1.9, compared to Pq (s) (dashed line) and Pff (s) 
(continuous line). B: parameter rj as a function of e/eF for = n/2. r s — 0.8 (Fermi glass, squares), 
r s — 6.3 (intermediate regime, full circles) and r s — 30 (Wigner crystal, diamonds). C (Ensemble 
average < \6\ >) and D (typical total longitudinal current /;) as a function of e/eF (same cases and 
symbols than in B). Inset of C: angular distribution P\g\ in the intermediate regime (r s — 6.3) for 
e/eF — 0.4 (dashed line), 1 (dotted line) and 1.4 — 1.9 (solid line), respectively. 



being expected in the thermodynamic limit. For the intermediate regime (r s = 6.3) one has 
quantum ergodicity and W-D statistics above a (low) excitation energy w cf, but the first 
excitations exhibit deviations from W-D statistics, indicating a break-down of a description 
in terms of weakly interacting quasi-particles and suggesting possible non chaotic collective 
excitations. 

The low energy excitations are characterized by a particular structure of the corresponding 
persistent currents when a flux <p is applied. The current /'") of the n th many-body level has 
a total longitudinal component given by 



ii n) (4>) = - 



dE n 






(3) 
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(n) 

The local current j^ flowing at the site i in the longitudinal direction is defined as 

j\f = 21m(* n \cl +14v c ix , iy exp(t0/X)|* n ), (4) 

where SSf n (n — 0,1,2,...) is the many-body wavefunction of energy E n . One defines the 
local transverse current j\ n t by a similar expression. The local current of the n th level 
flowing at the site i is a two-component vector j^ n \i) — (ii? >Jit) characterized by its angle 

Of 1 ' = arctan( j\ y / j\ n t ' ) and its absolute value j\ n ' = \j^(i)\- 

The inset of Fig. |l| C gives the distribution P| | of \8\ in the intermediate regime for 
increasing excitation energies (Pg is symmetric when 8 — > —0). For e/ep > 1.4 where W-D 
rigidity sets in, P\g\ is uniform, the flow is randomly oriented. For e/eF < 1.4 where deviations 
from W-D statistics can be seen, P\g\ is non uniform, the flow being mainly paramagnetic at 
e/eF ~ 0.4 and diamagnetic at e/eF « 1. The average of \9\ is shown in Fig.[j] C: values < 7r/2 or 
> 7r/2 correspond to paramagnetic or diamagnetic responses respectively. In the Fermi glass, 
the flow is weakly paramagnetic and becomes randomly oriented for e > ep. In the pinned 
Wigner crystal, a strongly paramagnetic response becomes randomly oriented at e > 3e^. 
The intermediate regime has a very specific behavior, the flow being mainly paramagnetic 
(e/eF < 0.7) then diamagnetic (e/eF ~ 1), before being randomly oriented at e/eF > 1-4 when 
quantum ergodicity sets in. The amplitudes of the longitudinal paramagnetic (probability c p ) 
and diamagnetic (probability Cd) total current Ii(E n w e) have log-normal distribution. The 
typical value /; = c p exp < ln|/ ; p | > — c^exp < ln|/ ; d | > is given in Fig. p] D. Below the 
ergodicity threshold of the intermediate regime, Ii is non random, being paramagnetic then 
diamagnetic in contrast to the Fermi and Wigner limits. The sign of /; can be explained 
[181 in the Wigner limit, where the </> dependent parts of the energies E n can be expanded in 
power of t/U . The leading contribution comes from the shortest Id paths enclosing (f> and the 
sign is given by Leggett's theorem (see Ref.||) for Id spinless fermions. The sign depends on 
the sample geometry and the parity of N when r s — > oo, independently of the microscopic 
disorder, but remains the same as far as r s > rf , becoming sample dependent in the Fermi 
glass only. 

Additional informations on the spectral statistics are given in Fig. [2| A (ij as a function 
of r s for levels of increasing excitation energies). There are two values of r s (sa 2 and « 30) 
where P(s) is invariant when e/ep < 3. Between those two values, the spectrum becomes 
more rigid as one increases e, while the opposite behavior is observed otherwise, providing 
another signature of the intermediate regime. Assuming a Gibbs-Boltzmann population of the 
many-body levels, we show in Fig. B the variation of 

«°w"ii>=( £ -j | fr?; ( \-,y ) w 

\ 2^ n ^p(-(E n - E )/kT) / 

as a function of r s at different temperatures. Another signature of the intermediate regime 
is obtained, the enhancement of the above expression disappearing as T — > Tp , while there is 
a weak increase in the Fermi and Wigner limit. We point out that the metallic (insulating) 
phase is identified in the transport measurements from a decrease (increase) of the conductivity 
when T increases and that the low temperature metallic behaviour disappears for T sa Tp- In 
the inset, one can see that enhancement and opposite temperature dependence disappear for 
a larger disorder (W/t — 15), indicating that the ratio W/t is another relevant parameter, a 
metallic behavior being not expected |19f for too strong disorders. 

The distribution Pg calculated for the ground state ( Fig. H A) shows how the local angles 
undergo a transition from a nearly uniform angular distribution (glass of currents randomly 
directed in the plane) towards a distribution strongly peaked in the longitudinal direction 
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Fig. 2. - W/t — 5 and <j> — 7r/2. A: Parameter jjasa function of r s for the first excitation so (circles), 
S2 — Si (squares) and Sio — S20 (diamonds). B: Typical longitudinal current for temperature in units 
of the Fermi temperature T/Tf = (circles), 0.4 (squares), 0.7 (diamonds), 1.1 (triangles up). Inset: 
the same at W/t — 15. 
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Fig. 3. - 4> — t/2. A: distribution Pg of the ground state local current angles for r a = (dotted 
line), r s — 6.3 (full line) and r a = 42 (dashed line). B: < \8\ > (open symbols, left scale) and 
amplitudes < In j > (filled symbols, right scale) as a function of r s for the ground state (circles) and 
for 1.4 < e/e F < 1.9 (triangles). 



(ordered flow of currents with a Id topology). Fig. B gives the average angle < \8\ > and 
the local current amplitude < lnj > when r s increases, both for the ground state and for 
e > cf'- The pattern of driven currents becomes ordered for r s values smaller than those at 
which j is suppressed due to charge crystallization. 

In summary, the local currents and the level statistics characterizing the low energy exci- 
tations confirm the existence of an intermediate regime in small clusters. The main issue is 
to know if this intermediate regime gives rise to a new quantum phase in the thermodynamic 
limit. A finite size scaling (FSS) study using exact diagonalization (larger L at the same 
density n e ) is beyond the ability of the most powerful available computers. One alternative 
could be to use variational trial wavefunctions justified for small and large r s and to compare 
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their relative energies using Monte Carlo method GO]. However, such approaches depend on 
the taken trial wavefunctions and can easily miss complex behaviors as the ones behind the 
beautiful FQHE physics, where the exact numerical solution |2^] of the three particle problem 
have been instructive. To study larger sizes, we have also developed approximate methods 
where the huge Hamiltonian is truncated in restricted sub-spaces using one body states pl| or 
Hartree-Fock orbitals |l^] . From those studies valid for small and large r s , but which become 
unreliable for intermediary values of r s , we know that the obtained precisions do not allow us 
to calculate very small quantities (as the local currents presented here). Numerical studies of 
those Coulomb systems are limited to exact diagonalizations of small clusters combined with 
approximate FSS approaches. For instance, a FSS study [Ell of the length £ characterizing the 
change of the ground state density after a local perturbation of the substrate has revealed the 
existence of a divergence when r s — > rf , giving the signature of a real quantum transition. In 
this case, we propose the local current j as a possible order parameter driving the transitions 
from the intermediate phase to both the Fermi glass (via its angle 9) and the pinned Wigner 
crystal (via its amplitude j). The melting of the Wigner crystal leads to an intermediate 
"hexatic" phase which has no translational order but still has orientational order (alignment 
of the local currents j). Shear modulus measurement |23] would be very useful, together with 
transport and compressibility measurements, to clarify these issues. We note that the disorder 
strength W should play a role which has not been fully explored (the study is mainly centered 
on kpl ~ 2.7). The role of the spin degrees of freedom will be discussed in details in a following 
study |2J] which does not lead us to fundamentally reconsider our conclusions. 
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